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ABSTRACT 

We investigate the clustering of dark energy within matter overdensities and voids. In particular, 
we derive an analytical expression for the dark energy density perturbations, which is valid both in 
the linear, quasi-linear and fully non-linear regime of structure formation. We also investigate the 
possibility of detecting such dark energy clustering through the ISW effect. In the case of uncoupled 
quintessence models, if the mass of the field is of order the Hubble scale today or smaller, dark 
energy fluctuations are always small compared to the matter density contrast. Even when the matter 
perturbations enter the non-linear regime, the dark energy perturbations remain linear. We find that 
virialised clusters and voids correspond to local overdensities in dark energy, with 5^/ (l+w) ^ 0(10^^) 
for voids, (50/(1 -I- w) ~ 0(10"'') for super-voids and ^^/(l + w) ^ O(10~^) for a typical virialised 
cluster. If voids with radii of 100 — 300 Mpc exist within the visible Universe then (5^ may be as large 
as 10~^(1 -|- w). Linear overdensities of matter and super-clusters generally correspond to local voids 
in dark energy; for a typical super-cluster: (5^/(1 + w) ^ 0{—lQ~^). The approach taken in this work 
could be straightforwardly extended to study the clustering of more general dark energy models. 
Subject headings: Cosmology: Theory, miscellaneous. Relativity. Galaxies: general, Large scale 
structure of the universe. 



1. INTRODUCTION 

It has been almost a decade since observations of 
Type la supernovae (SNe la) were first found to sup- 
port the notion that our universe is c urrently undergo- 
ing a phase of accelera ted expansion (jRiess et al.l 119981 : 
iPerlmutter et all 119990 . Since that time, evidence in 
favour of this accelerated expansion has strengthened 
significantl y as the result of further SNe la ob servations 
(jRiess et a h 2004, 2006a^b; Wood-V asev et al.l i2007). im- 
proved measure ments of the cosm i c mic rowave back- 
ground (CMB) (jSpergel et a l ] l2003l [200l and surveys 
of lar ge scale structure (LS S) (jAdelman-McCarthv et al.l 
I2006t iTegmark et al.l [20060 . The precise cause of this 
late-time acceleration, however, remains unknown. 

If general relativity is accurate on astrophysical scales, 
then the assumptions of large-scale homogeneity and 
isotropy require that the agent responsible for the uni- 
verse's acceleration, dubbed 'dark energy', behave cos- 
mologically as a fluid with negative pressure. The 
standard model of particle physics predicts only one 
such fluid: the vacuum energy, or cosmological con- 
stant, for which the pressure, p, is always equal to mi- 
nus the energy density, p. However, if the vacuum en- 
ergy density is indeed non-zero then it is generally ex- 
pected to be of the order of Mpj, where Mp\ is the 
Planck-mass. This is some 120 orders of magnitude 
larger than the observed dark energy density. A num- 
ber of proposals have therefore been made in the lit- 
erature for models in which dark energy is dynamical 
and associated with some new form of energy (W etterichl 
119881 iPeebles fc Ratral 119880 . In these models, the 
present small size of the effective cosmological constant 
is by-product of the age of Universe. The accelera- 
tion of the Universe might alternatively be explained 
by modifying General Relativity rather than postulating 



a new fo rm of energy (s ee iDvali, Gabadadze fc Porratil 
'2000; Ama rzguioui et al . 2006; Noiiri & Odintsov 200l 
j Capozziello. Cardone fc Troisi ,2005: ; Koivisto fc Motal 
boOTallbl ). It also been suggested that the Universe is 
not accelerating at all, but that a large local inhomo- 
geneity prevents the SNe la data from being correctly 
interpreted in ter ms of a homogeneous and isotropic cos- 
molo gical model (|Kolb Matarrese fc Riottoll2006l : iMoffalj 
I2OO60 . This said, dynamical dark energy (hereafter 
DDE) is by far the most popular candidate to explain the 
current astronomical data. In the simplest DDE models, 
known as quintessence, dark energy is associated with the 
energy density of a scalar field with a canonical kinetic 
structure. 

If dark energy does indeed exist then astronomical ob- 
servations presently provide us with only hints as to its 
nature. We know that today it represents about 70% 
of the total energy density of the Universe, and that its 
equation of state (EoS) parameter, w = p /p, is fairly 
close to —1: w = — 1 ± 0.1 for z < 1, (jRiess et al.l 
[200 6a). Furthermore, if matter dominates the expansion 
of the Universe for z > 1.8 then w{z > 1) = -0.8tQi 
(jRiess et al.l l2006aO . Hence, dark energy has negative 
pressure at higher redshifts with a 98% confidence. These 
bounds on w are entirely consistent with a pure cosmo- 
logical constant. Whilst detecting either w ^ —1 or 
dw/dz ^ would rule out a cosmological constant, in 
many DDE models significant deviations from w = — 1 
only occur at early times and as such would be difficult to 
detect. The late time, background cosmology predicted 
by many DDE models is therefore very similar to that of 
a universe with a true cosmological constant. 

DDE models generally cease to mimic a cosmo- 
logical constant in inhomogeneous backgrounds or 
when one considers cosmological perturbation the- 
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ory. In particular, a number of authors have stud- 
ied the ef fect of DDE on the formation of large scale 
structure (ICorasaniti. Giannanto nio fc Melchiorri' '2005'; 
Hannestad fc Mortselll 120021: iBrookficld ct al. 2006b; 



Doran. Robbers fc WetterichI 120071) . In the vast major 



ity of these works the energy density of dark energy is 
taken to be homogeneous i.e. it is assumed that DDE 
does not cluster. The extent to which this assump- 
tion of homogeneity is valid has been the subject of 
some interest and much debate in the literature. In 
an inhomogeneous background, the DDE energy den- 
sity and EoS parameter should exhibit some spatial vari- 
ations. The key issue of how large these variations 
should be is however far from settled, pa rticularly when 
the matte r perturb ation goes non-linear (Maor & Lahav 
2005; Mota fc van dc Br uck 2004: Linder fc Wh ite 2005; 
Bartclmann, Doran fc WetterichI 120051: iNunes fc Motel 
2006: .Abramo et all I2007D . In this paper we attempt 



to settle this issue for uncoupled quintessence models by 
deriving an analytical expression for the dark energy con- 
trast, (50, in by presence of a matter perturbation. Im- 
portantly, we do not constrain the matter perturbation 
to be small (i.e. linear). 

Clustering of DDE over scales smaller 
than about 100 Mpc has been th e sub- 
a number of recent articles (iCaimmil 
iBalaguera-Antolmez. Mota fc Nowakowskij 



of 



[M 



ammi 



20051: iPercivall l2005r 



iNunes. da Silva fc Aghanimf 



IWan 
TgO^ 

Balaguera-Antolmez. Mota fc Nowakowskil l2007t ). 
Most attention has been focussed on models 
which the DDE couples to baryonic and/or 
matter since DDE clusterin g is expec ted to 
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Brookfield et al.' '2006a'; 'Manera fc Motal 12006 
Pettorino, Baccigalupi fc Mangano 2009). In other 
works, a more phenomenological approach has been 
taken and the observables associated with DDE cluster- 
ing have been parametrized. In particular, it has been 
shown that inhomogeneities in dark energy could pro- 
duce detectable signatu res on the CMB (WcU er fc Lewis! 
[200l iKoivisto fc M ota 200l). 

Recently. iDutta fc Maor (|2007t l considered the growth 
of spatial DDE perturbations in the absence of any mat- 
ter coupling. They considered only those circumstances 
where both the density contrast of matter, Sm, and that 
of the DDE, S^, were small enough to be treated as a lin- 
earized perturbations about a homogeneous and isotropic 
cosmological background. They studied the simplest 
class of quintessence models where dark energy is associ- 
ated with the slow-roll of a scalar field (j> down a poten- 
tial V {((>); 4> is minimally coupled to gravity. The authors 
linearized the full field equations for 0, the matter and 
the metric and solved the resulting system numerically. 
Intriguingly they found that, at late times, a local over- 
density of matter corresponded to a local under-density, 
or void, of dark energy. Conversely, a void in the mat- 
ter was seen to produce a local DDE overdcnsity. Al- 
though the DDE density contrast, (5^, is initially very 
small compared to the matter density contrast, Sm, they 
found that, when Sm ~ \S^\ 0(10"^). S^ was 

also observed to be growing more quickly than Sm at late 
times. Their results suggest that DDE clustering may in- 
duce a positive correction to the value of (l-f w) of more 



than 10% at the centre of an inhomogeneity with prop- 
erties similar to that of the local supercluster (hereafter 
LSC). If accurate, the results of iDutta fc MaoA (|2007D 
imply that any deviations from w = —1 would be signifi- 
cantly amplified by the presence of a local overdensity of 
matter. Furthermore, they suggest that DDE clustering 
might be relatively strong when the matter perturbation 
g oes non-linear. 

IDutta fc Maori ()2007l ) studied DDE clustering when 
the matter perturbation is growing in the linear regime, 
which is only accurate when Sm ^ 1. They considered 
the evolution S^ for a cluster of matter with initial den- 
sity profile (at z = 35) of Sm = Aexp{—r^/a^) where 
A = 0.1 and cr 0.01H~'^ with Hi being the initial 
value of the Hubble parameter. However, in the absence 
of any DDE, the linear approximation would generally 
cease to be accurate when z « 2.8; additionally the per- 
turbation would be expected to turnaround when z « 1.0 
and virialise when z ~ 0.3. It is therefore far from clear 
wh ether or not th e shar p late-time growth in S^ found 
bv IDutta fc Maori (|2007t l is indeed a physical effect, or 
just a result of using linearized field equations outside of 
realm in which they are valid. Although one might well 
expect there exist some mapping between the linear and 
non-linear regimes as it happens in an Einstein-de Sitter 
Universe. 

In this paper w e investigate a simil ar problem to 
that considered bv IDutta fc Maori (|2007D . i.e. the clus- 
tering of uncoupled quintessence on sub-horizon scales. 
There are, however, two important differences between 
our approach t o thi s problem and the one taken by 
IDutta fc Maori (|2007f ): 

1. Firstly, we do not use numerical simulations but 
instead use the method of matched asymptotic ex- 
pansions (MAEs) to develop a analytical approxi- 
mation to (50. 

2. Secondly, we do not require the density contrast of 
the matter perturbation, Smi to be small. Indeed, 
our analysis and results remain valid even after the 
virialisation of the matter overdensity. 

Our approximation for S^ is accurate provided that Sm is 
only non-linear (> 0(1)) on sub-horizon scales. This pro- 
vision is consistent with observations. For simplicity we 
take the matter perturbation to be spherical symmetric. 
We also require that gravity is suitably weak whenever 
the inhomogeneity is non-linear i.e. GM/R <C 1. We re- 
state these requirements in a rigorous fashion later, how- 
ever they are essentially equivalent to the statement that 
gravity is approximately Newtonian over scales smaller 
thani?"^. Since DDE clustering in the linear reg i me ha s 
been dealt with in great detail bv IDutta fc Maoil (|2007D . 
our main focus in this paper is on what occurs when the 
matter perturbation goes non-linear. 

This paper is organized as follows: in Section [2] we de- 
scribe our model for an inhomogeneous spacetime and 
the DDE. We state the equations that must be satisfied 
by the metric quantities and the dark energy scalar field 
(f). In Section [3] we introduce the method of matched 
asymptotic expansions (MAEs). This method relies the 
existence of locally small parameters, and we state what 
these are for our model and interpret them physically. 
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We also note what constraints the smallness of these pa- 
rameters places upon our analysis. In Section |4] we ap- 
ply the method of MAEs to the evolution of dynamical 
dark energy perturbations. We derive a simple equa- 
tion for the DDE density contrast, S^, in terms of the 
peculiar velocity of matter particles, Sv^ and the mat- 
ter density contrast, 5m- Importantly this equation is 
equally as valid for Sm > 0(1) ('the non-linear regime') 
as it is when Sm <C 0(1) ('the linear regime'). In Sec- 
tion [5] we use our results to study the evolution of 5^ in 
the hnear {Sm ^ 1), quasi-hnear {Sm ~ 0(1) — 0(10)) 
and fully non-linear {6m ^ 1) regimes. We compare our 
analytical result s with those found n umerically in the 
hnear regime bv lDutta fc Maoil (|2007( ). In Section [HI we 
consider the spatial profile of the DDE density contrast 
for realistic astrophysical inhomogcncities such as voids, 
supervoids, clusters and superclusters. We conclude in 
Section [7] with a discussion of our results and their ob- 
servational implications. We also note how our analysis 
might be extended to include even more general DDE 
models. 

Throughout the paper we use units where c = 1 . 

2. THE MODEL 

2.1. Geometrical Set-Up 

Our aim in this paper is to derive an expression for 

the DDE density contrast, = Se^'^'' /ei'^\ inside an 
over-density, or under-density, of matter. Henceforth we 
use e*-*^ to represent the energy density of a component 
i. We are particularly interested in those cases where 

the matter perturbation is non- linear i.e. fe^^Vec™^ > 
0(1). Although we aim to remain suitably general in our 
treatment of the density perturbation, we do make the 
following simplifying assumptions: 

• We assume spherical symmetry. We briefly discuss 
the extent to which the relaxation of this assump- 
tion would affect our results in Sect ion [7] below, and 
conclude that the qualitative nature of our findings 
would be unaffected. 

• We define a 'physical radial coordinate', R, by the 
requirement that a spherical surface with physical 
radius R has surface area AnR^. We assume that 
all curvature invariants are regular at i? = i.e. we 
do not consider those cases in which there is a cen- 
tral black hole. We argue below, however, that our 
results are still accurate even when there is a cen- 
tral black-hole, provided they are only applied at 
radii that are large compared to the Schwarzschild 
radius of the black hole. 

• Finally we assume that for radii smaller than some 
Rq ^ gravity is suitably weak in the inhomo- 
geneous region; H is the Hubble parameter of the 
background spacetime. We define what we mean 
by weak rigorously in Section [3l For radii R > Rq, 
we require that the matter density contrast is < 1. 
This assumption holds for most realistic models of 
collapsing overdensities provided that the radius 
the overdense region is less than about 0.1/ H. 

2.2. Einstein's Equations 



We take the matter content of the Universe to be a mix 
of irrotational dust and dynamical dark energy, which 
is described by a scalar field </). For simplicity we re- 
strict ourselves to considering only spherically symmet- 
ric spacetimes for which the most general line element in 
comoving coordinates is: 

ds^ = dt^ - U{t, r) dr^ - R{t, ( d9^ + sin^ 6 dif^) . 
We make the definitions U{t,r) = R'^j.{t,r)/Q{t,r) and 
k{t,r) = 1 — Q{r,t). This coordinate choice is unique 
np to t ^ t + to, and r r'{r). With these definitions 
the 2-sphcres {t, r} = const have surface area AuR^it, r), 
and in this sense R{t, r) represents the 'physical radial 
coordinate'. The energy- momentum tensor of pressure- 
less dust is given by 



(m) 



ab 



dia; 



lag (e(™\ 0,0,0 



and the energy-momentum tensor of the scalar field, 
is 



(0) 



ab 



da(, 



'Jb(f> ~ gab ( 



=0 - y(0) 



The Einstein equations for this metric read: 



Gab — Rab 



2^9 ab 



K T 



T, 



(1) 



where k = 8tt. The tt-component of Eq. ([T|) gives: 

{R^t.R + kR) = R,rR^Ke^"'^ + R.rR^n ( i./.^ (2) 



and the ir-component of the Einstein equations is: 

^^ = -K0(/',r- (3) 

The rr-component of Eq. ([T]) reads: 

(R^tR + kR) +Q^tR^ (4) 



+ - ] R,R^ 



^From r(™)'';,.a = 0, it follows that: 

where F(r) is an arbitrary constant of integration. If we 
define r so that R{r, ti) = r for some t = ti, then F{r) = 
2M^^r (r) I \/Q{r, ti) where Mi (r) is the mass inside a shell 
of radius r at t — ti. 

iFrom T^'>''>''b;a = one obtains: 

- 00 = V^{^). (6) 
Eq. © is subject to the boundary conditions: i/i^r = 
at J? = 0, and 4>c{t) — hm^-foo 0('', i), where 0c(i) is the 
solution of Eq. ^ in the cosmological background. We 
make the following definitions: 



2M{r,t)^ F{r)^Q{r,t)dr, 

5e^''^-\^' + i^fr + V{4>)-e,{t), 
5eM=Se^'^^-Rj^, 

R r 
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where ro(t) is defined by R{ro{t),t) = 0. We also define: 

Integrating Eqs. ([3]) and ^ and using Eqs. (O and © 
we find tfiat 



K 



, 2M (r, t) 1 , / . ,^ , , 



+ 



R 



ro{t) 



R 3 

R^r{z,t)R^{z,t)5e^'''\z,t) dz, (7) 



and 



M(r, t) 1 



i?2 



6 



(8) 



Jro(t) ' 2 

As i? ^ oo we must recover tlie FRW background cos- 
mology which implies: 

lim A^(r,^)^l«:4™)(^o)a'(^o)r^ 

r — *oo Z 

lim fc(r, t) ^ kor'^ , 

r — ^oo 

where to is an arbitrary time, a{t) is the scale factor 
of the FRW background, and ko is the curvature of the 
background. In line with current observations, and be- 
cause it greatly simplifies the calculations, we take the 
cosmological background to be flat and set ko — 0. We 
do not attempt to solve the Einstein or DDE equations 
exactly, but instead we develop asymptotic approxima- 
tions to the true solutions which are accurate so long as 
a number of parameters, which we define and interpret 
below, remain small. 

3. SMALL PARAMETERS AND THE MATCHED 
ASYMPTOTIC EXPANSIONS. 

One can think of the small parameters approach and 
the matched asymptotic expansions as an expansion in 
the Newtonian potential. Although it is perhaps fairer 
to say that is closest to the context of General Relativ- 
ity. The differences are not obvious at leading order but 
at next to leading order one will see that what appears 
in the equation for the acceleration is not necessarily 
what one would expect to appear in the Non-relativistic 
regime. 

As stated above in Section [2?T1 we require that gravity 
is 'suitably weak' inside the density perturbation. By 
suitably weak we mean that 



where 



\Si\ « 1, 
Si=6v = R^t - HR. 



Physically R^t is the velocity of a particle of dust in {t, R} 
coordinates, where R is the physical radial coordinate, 
and HR is the velocity that such a particle would have in 
the cosmological background. The small parameter Si = 
Sv is therefore the peculiar velocity of a matter particle 
with respect to the cosmological background. Recall that 
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Fig. 1. — Sketch showing the interior, exterior and intermedi- 
ate regions. The interior is defined to be the region where the 
inner approximation is valid, and the exterior is the region where 
the outer approximation is applicable. In order for the matching 
procedure to work, we must require that the exterior and interior 
regions overlap in some intermediate region. We must also require 
the exterior and interior regions are contiguous, so that one may 
move smoothly from the exterior, through the intermediate region, 
to the interior. 



we use units where c = 1, hence the first expression above 
is in fact the dimensionless quantity Sv/c <^ 1. Our 
assumption requires that whenever HR > 0{1), we have 
\di/HR\ <C 1. This condition is equivalent to requiring 
that the mean matter density contrast, 6{R,t), inside 
the sphere with radius R ~ 0{1/H) or greater, is small 
enough that it may be treated as a linear perturbation 
about the cosmological background. The mean matter 
perturbation is however allowed to be non-linear, 6 > 
0(1), when HR < 1 provided |i?,t| < 1. We also define 



(m) 



1 



) = 



We require ^2 ^ 1- This again require that S„i > 1 
only when HR <C 1 i.e. the matter perturbation is only 
non-linear on sub-horizon scales. We note that 62 ^ 
0{{R'^t + k)- H'^R^), and so generally |<5i | < 1 provided 
52 < 1. 

In what follows we assume that i5i and 82 are small 
everywhere. These assumptions can be checked once one 
specifies initial conditions for the matter overdensity, but 
they generally hold very well whenever the scale of the 
inhomogeneous region is < 0.1/ H. 

We now define two over-lapping regions which wc shall 
refer to as the interior and the exterior. 

3.1. The Interior Region 
The interior region is defined by 
S3 = HR<^ 1. 

We define the inner limit of a quantity F{r,t) to be 
Lint{F) = Vnns^^o F(r^t). The inner limit can be imag- 
ined as the limit in which the cosmological background 
density of matter is taken to zero, and the cosmologi- 
cal horizon is taken to infinity. In the interior we con- 
struct asymptotic approximations to quantities in the 
limit HR 0. We also expand in (5i and 82. We re- 
fer to this as the inner approximation. 
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3.2. The Exterior Region 

In the exterior region, spacetime is required to be 
homogeneous and isotropic at leading order. In other 
words, the matter perturbation must be hnear in the ex- 
terior: 

^4 = \Si/63\ < 1,(55 = \S,n\ < 1. 

From the form of R^t, it is clear that 85 ^ implies 
64 0. We define the exterior limit of a quantity F{r, t) 
to be Lext{F) = limg^^o F{r,t). As in the interior limit 
we also expand in i5i and 62- We construct asymptotic 
approximations to quantities in the exterior region in this 
exterior limit and refer to them as the outer approxima- 
tions. 

3.3. Matching and the Intermediate Region 

We are primarily concerned with the behaviour of the 
DDE in the interior region. In Section [H we find the 
inner approximation to by solving Eq. ([6]) order by 
order in the interior limit. We cannot, however, ap- 
ply both of the boundary conditions on directly to 
the inner approximation. This is because the condition 
lim/j^oo (f' = 4'cit) must be applied at i? = 00, a point 
that is very clearly in the exterior and not in the in- 
terior region. As a result, the inner approximation will 
contain ambiguous constants of integration. Fortunately, 
this ambiguity can be lifted by matching the inner and 
outer approximations to (f> if there exists some intermedi- 
ate region where both approximations are simultaneously 
valid. This matching of the inner approximation to the 
outer one is referred to as the method of matched asymp- 
totic expansions (MAEs). It relies on the fact that, in 
any given region, the asymptotic expansion of a quantity 
is unique (for a proof see .Hin ch 1991). Thus if the inner 
and outer approximations of (j) are both valid in the in- 
termediate region, they must be equal in that region. For 
the method of MAE to be applicable we must, of course, 
require that an intermediate region exists. A necessary 
condition for an intermediate region to exist is that for 
some range of R: 

HR < 1, \-i52/ne^r\t)R^\ = |3,5™(t)/K| 1. 

This becomes a necessary and sufficient condition if the 
only boundary of the interior region is in the exterior 
one, and vice versa (see Figure [T] for an illustration). If 
the scale of the inhomogeneity, i?o, is taken to be the 
largest value of R for which |3i52/Kei'"^(t)i?^| > 0.3, and 
then then an intermediate region generally exists pro- 
vided that Rq < O.IH. 

In resume, in this section we have considered spher- 
ically symmetric backgrounds and far from black hole 
horizons. We have basically performed an expansion 
in the Newtonian potential at leading order, at least 
as far as the evolution of scalar field perturbations are 
concerned. Notice, however, that this method can be 
straightforwardly ex tended to allow for dey iations from 
spherical symmetry (jShaw fc Barrowir2006al lbll3) . 

To be clear, we should point out that the expansion 
in Si and 62, which is a lot like expanding in the New- 
tonian potential, is not really the key to our method. 
That is just a simplification. The key is that we take: 
S3 — HR <C 1 in the interior region and 5^ ^ 1 in 
the exterior region. That is, everything is linear in those 



two regions. We then assume that there is an interme- 
diate region where both of these conditions hold. Which 
implies that in the intermediate region: \Si\ <C 1 and 
1(^2 1 ^ 1. So the assumption that Si, 82 <C 1 everywhere 
ensures that we have an intermediate region. We could 
in fact relax this but the calculations would become more 
difficult: often still doable, but in almost all physically 
interesting cases totally unnecessary. For instance the 
relaxation to \Si\ < 1 and \S2\ ^ 1, is very straightfor- 
ward and allow s one to g o all the way up to a black hole 
horizon ( Shaw fc Barrow. 2006 a. bx:). 

4. EVOLUTION OF DYNAMICAL DARK ENERGY 
PERTURBATIONS 

In this section we find an asymptotic approximation 
to the DDE density contrast, S^. The discussion is fairly 
technical, and readers more interested in the results than 
the machinery used to derive them may prefer to focus on 
the statement, discussion and application of our results 
in Section 14.31 and following. 

The DDE is described by the field <j) which satisfies: 

We solve this equation by constructing an asymptotic 
approximation to (f) in the small parameter Si. We write: 

(j)^(j)cit) + Scj){t,r){l+OiSi)), 

where (50 ^ 0{Si). Before solving for (j>, we make the 
dependence of the metric on Si explicit by transforming 
to a new radial coordinate, p — R{t,r)/a(t), where a(t) 
is the scale factor of the Friedmann-Robertson- Walker 
(FRW) cosmological background. In {t, p} coordinates 
the metric is: 



ds 



2 _ {l-k{r,t)~Sj{r,t)) dt'^ Ji{r, t)a{t) At dp 
1 - k{r, t) ^ 1 - k{r, t) 



a^{t) dp^ 



-a'^p'^ {d0^ +sin^e dip} . 



1 - k{r, t) 
Consider — □(^c(^) in this metric: 



(9) 



1 u^ Mt) d I 2 



-□0c(^)=K.(0.) + ^t'^^(^)-«(t)p^ap 



p'Si 



where Q is evaluated at constant r as opposed to constant 
p or R and is thus given by Eq. ([3]). It follows that S(j) 
satisfies: 



-US^{\ + 0{Si)) = Cl/^(<^) - 1/0(<^c)) (10) 

We now solve for 5(j) in both the interior and exterior 
regions. 



4.1. Inner Approximation 

We note that \k^IR^ < 0((5|), and that \k\ - 
0{S3,Sf,S2), so that 

- ^dt {aHt)S^,t) + ip'Scjy^p) ^ {V^cj,) - V^^,)) a'p' 



+Mt)a (p'(5i) + OiSiS<j), S2S<j), S 
In many cases S(j) is small enough so that 

iV^icj)) - K0('^c))aV' - V^^i(bc)R^S<j), 
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and Vcficf, ^ 0{H^) so that (j>c{t) evolves over cosmological 
thiic-scalcs. In these cases it is clear that: 

and so this term can be dropped at leading order. Al- 
though it will often be the case that T^</,(0) — T^0(0c) ~ 
y,<ii<ii{4>c)5(l) we do not have to require such a strong as- 
sumption in order to justify ignoring the effect of the 
potential to leading order in the interior approximation. 
Instead, we make the far weaker assumption that: 

\V^{^c + S^) - K0(0c)|i?V'5'/' ~ 0{6i,d2, S3). (11) 

A physically-viable DDE theory for which this require- 
ment does not hold would be difficult to construct. For 
instance, if this condition did not hold then even if, at one 
instant, ^^(^c) ~ 0{H^) or smaller, as must be required 
for 4>c to evolve over cosmological time scales, shortly 
afterwards small changes in 0c would cause V^^{4'c) to 
grow to be many orders of magnitude greater than . 
This would result in (pc evolving over time-scales that are 
much shorter than the Hubble time. The only DDE the- 
ories that we can imagine that could accommodate this 
and still be compatible with observations would involve 
(j) being held almost completely fixed at a minimum of 
V{(t)c), in which case the DDE would be almost entirely 
indistinguishable from a cosmological constant. 

To leading order in the interior, we would actually be 
justified in setting a(t) = const as dp = HR = <C 
1. However, since it might not be completely clear that 
these terms may be ignored, and because we can solve the 
5(j) equation without having to ignore them, we continue 
to include them. 

We define d?7 = dt/a{t) and u — aScj), and note that 
p^a^^^/a = HR^ + H^R^ - ©(Jg), so that: 

+ ^ (p'u,p) )j (1 + 0(<5i, <52, <^3)) ~ 



;j,,[t)^[a{t)pH,)^^. 
which has as a solution 
<50(t,p)' 
^{■n ~ p) _|_ + p) 

a{t)p 0'{i)p 



(12) 



2 Jo 7-1 ait)X{p,p',s) 



+ 0{6i6cj),S25(b,d35(b), (13) 



where /(r;,p) = {p^,(t,p)) and X(p,p',s) = 



Vp^ 



p'"^ — 2pp's] T and Q represent waves in the scalar 



field and we discuss them further momentarily. To go 
further, we make the following reasonable assumptions 
about the behaviour of b\. We assume that there exists 
some i?_2, in the interior region, such that for R > i?_2 
in the interior region, 5i decreases faster than 1/R^ as 
R 00; / then decreases faster than 1/R'^ for R > i?_2- 
This assumption ensures that dominant contributions to 
the integral in Eq. comes from values of p' that 

are well inside the interior region. Note that generally 
f{ri,p) varies over conformal times-scales of order 1/aH 
and Si/aSi. Provided that Si is not momentarily zero, 
these scales are much larger than X < V2-R-2/a so we 
can Taylor expand /(?] — X, p') as 

/(77 - X, p') - /(77, p') + /,,(,7, p')X + 0(/,,,X^), 



where the neglected terms are 0(5^ , (5|, Jif^s, (^2) times 
f{ri,p'). After some algebra and integration wc then ar- 
rive at: 

6<Pit, R)il + (6^,62, S3)) ~ Mt) C ^^i*^ R') d^' 

I y^iv-p)+G{v+p) (|24) 

where the previous assumption relating to the Taylor ex- 
pansion of / can be seen to be equivalent to the statement 
that Si{t,R)R'^ <C Si. This assumption will generally 
break down if there is some initial instant, t = ti say, 
when (Si = 0. M t ti however, jii?^ <C Sx is equiva- 
lent to S2 ^ 1, which is certainly the case in the 
interior region. 

Generally, the functions T and Q are related to the 
initial conditions on Sep. Because we have required 
R-2 ^ our analysis will break down at early times 
since generally HR-2 — > 00 as t — > 0. As a result, we 
cannot generally determine T and Q by simply speci- 
fying some initial conditions on(50, at t = ti say, and 
applying them to the interior approximation since the 
interior approximation by not be valid when t = ti. The 
requirement that 6(f>{R = 0) can however be applied to 
the interior approximation to give J^{r]) — —Q{r]). To 
the order at which we work, J-'{ri) itself should be fixed 
by matching to the outer approximation. 

4.2. Exterior Solution and Matching 

In order for cpc = linifl,_^oo '/'(i, R) we must have 
S(j){t^ i?) — > as -R — > 00. In the exterior region, space- 
time is FRW to leading order in and 82, and the 
equation describing for the leading order behaviour of 
u = a{t)S4> reads: 



^{p\,)^{V4 



-H- H'^] a^u + 



1 



K^'^{t)a'' pu^p 



Mt)a^ 



(p'Si)^ 



(15) 



2 - ■ p^ 

where dry — At/a{t). Solving this equation is far from 
straightforward. Fortunately, however, our main inter- 
est in not in how 5(j) behaves in the exterior region but in 
interior region where the perturbation in the matter den- 
sity can be non-linear, and so it is not necessary to solve 
Eq. (|15p in the exterior region so long as we know how 
its solutions behave in some intermediate region where 
S3 = HR ^ 1. We require that 50 — > as the mat- 
ter perturbation is removed i.e. Si — > 0. We have also 
required that drop off faster than 1/R for R > R_2 
where i?_2 is in the interior region i.e. HR^2 ^1- In 
this intermediate region, the solutions of Eq. ([15]) for 
which 50 — > as i? ^ 00 have the following form: 

50(i,i?)-0c(O / Si{t,R')dR' 

O((5i50,53(50,52<50), (16) 
where C{rj) ^ 0{1/H) or greater. We have assumed that 
t is larger compared to any initial instant when one one 
or more of Si — 0, S(p = or Scj) = hold. With these 
assumptions, matching to the interior region gives: 

Kirj)^J^{7^)+g{r^)^0, 

and 

^^Av) 



Kjy) 
R 



a{rj) 



- ~Mt) / Siit,R')dR'. 
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We have required drop off of 6i raeans that, to leading or- 
der, we can set C — oo. This impHes that the 2!F^r]/a{r]) 
is sub-leading order in the interior approximation and so 
may be neglected. 
Thus to leading order in the interior region we find: 

6(l>{t,R)^Mt) / Sv{R\t)dR'{l + 0{Su 62,63)). 



Where Sv — R^t — HR is the peculiar radial velocity 
of the matter particles relative to the expansion of the 
background Universe. This approximation is valid in the 
interior region at late times compared to any initial in- 
stance when 5(j) = Q and/or 50 = which is equivalent 
to \5v\R^ < \5v\R < \5v\. Note that, to this order, the 
surfaces of constant are surfaces of constant: 

/•OC 

t^^t- / 6vdR'. 
Jr 

We use this asymptotic approximation for Scf) to evaluate 
the DDE density contrast on surfaces of constant t in 
the matter rest frame. In the rest frame of the matter 
particles: 

sW^^-^^ + ^c^% + V{<p), 
and the local inhomogeneity in e^"^) is therefore: 
fc(«^02(^W_3^ f Sv{R',t)dR' + {Sv{R,t))^ 



j" [Ru -iH + H^)R'^ di?'^ 



(17) 



The corrections terms are 0{6i, 62,53) times smaller than 
the leading order term. Since 0^ = (1 + ii;)£c'^^ we have 

S^^{l + w)(-3H I dv{R\t)dR' + {Sv{R,t))^ 

+ J (^R,tt-{H + H^)R'^ dR'^, 

R,tt is the acceleration of the matter particles in {t, R) 
coordinates. Eq. ^ gives: 



It It a — ll it — II It 



5M 



R 

where 5M{R,t) = M{R,t) - \VL„,H'^R^ = \n„,H^RH 
is the mass contrast at time t inside the sphere with 
physical radius R. 

4.3. Discussion 

In this section we found that the DDE density contrast 
can be expressed, to leading order, in terms of the radial 
peculiar velocity, 5v , the mean density contrast, 5 and 
the unperturbed DDE equation of state parameter w: 

5^ ^ {1 + w) (^3H J 5v{R',t)dR' + (Sv)^ 

-ln„,H^ f 5{R',t)R'dR') . (18) 



This expression for 5^ is our main result and it is valid 
for both linear and non-linear sub-horizon matter per- 
turbations. In order to evaluate the above expression, 
one must simply specify the peculiar velocity, 5v,_of the 
matter particles and the matter density contrast 5. This 
exterior is valid in the interior region at late times com- 
pared to any initial instance when 5(1) = and/or 5(f> = 
i.e. \5v\R^ < \5v\R^ \5v\. 

Note that H and w will depend on the mass of the field 
m, so in that sense, equation (llSp is not completely model 
independent. However, all the model dependent terms 
come from the homogeneous and isotropic background 
cosmology. So the expression for 5^ is model independent 
provided mRc <C 1. Which is true if to ~ 0{H) and 
HRc <C 1 as assumed. 

The key point is that large deviations in 5^ from the 
prediction using a LCDM background could only occur if 
5v and 5m were predicted to change by an order of mag- 
nitude or more. As far as we know, however, structure 
formation is compatible with LCDM over the scales of 
clusters and so it seems unlikely that such large devia- 
tions occur. 

As we shall see below the DDE perturbation is gen- 
erally small compared to the matter perturbation, so to 
leading order one may neglect any DDE perturbations 
when evaluating 5 and 5v. 

It is interesting to note that the sign of 5^ depends 
on the relative magnitudes of 5 and 5v: From equation 
(fT8|) one deduces that the formation of local overdensities 
or local voids of dark energy arise from two competing 
effects: 

• The Drag effect: associated with deviations from 
5v/HR — 0. If one has an overdense region, 
then it expands mor e slowly. In agreement with 
iDutta fc Mao^ (|2007D . we find this drag effect to 
produce a local underdensity of dark energy. 

• The Pull effect: associated with deviations from 
i5m = 0. An overdensity of matter pulls the dark 
energy (and everything else) towards it. So matter 
and fields tend to clump around it. As a result, 
the pull effect means that an overdensity of matter 
creates an overdensity of DE. 

So there is one effect that pulls the dark energy in and 
another that pushes it out. As we will see, in the linear 
regime the drag effect grows more quickly than the pull 
effect, so one gets a dark energy void. The onset of non- 
linear structure formation, however, causes \5v/HR\ to 
reach a maximum and then to tend to 1 at late times, 
whereas 5m just keeps on growing. So the pull effect dom- 
inates at late times, so one ends up with an overdensity 
of dark energy. 

We have assumed above that the mass of the DDE 
scalar field is small compared to the inverse length scale 
of the cluster. Provided this is 5^ in uncoupled DDE 
theories is independent of the details of theory describing 
the dark energy. 

In the next section we use Eq. pB]) to evaluate 5^ 
in both the linear and non-linear regimes. We find that 
\5^[R = 0)1 ~ 0{VlmH'^Rl5ciust), where Rc is the radial 
scale of the cluster, and 5ciust is the mean matter density 
contrast va R < Rn. 
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5. CLUSTERING OF DYNAMICAL DARK ENERGY 

5.1. Dynamical Dark Energy Clustering in the Linear 

Regime 

Although one of our main aim in this work was to study 
DDE clustering in the non-linear regime, our analysis 
is also perfectly valid in the linear regime, i.e. when 
(5m 1, provided that the matter inhomogeneity is sub- 
horizon at the instant when we wish to evaluate 6^. We 
found above that is given by Eq. ((TS]). In the linear 
regime \Sv\ ^ HR, and so the term proportional Sv'^ 
in Eq. ^TE\i is small compared to the other two terms 
and should be neglected. 5m{r,t) is the matter density 
contrast, and we define Sk{t) to be its Fourier transform. 
In the linear regime (in comoving-coordinates): 



Evolution of central density contrasts In ttie linear regime 



and so: 



Svk = - 



-V • (5v, 
iHf{a)aSk 



where /(a) = dlni5fc/dlna and Svk is the Fourier trans- 
form of Sv. We recognize that in,nHHR/2 = 5M/R^ = 
dS^/dR where 5$ is the leading order perturbation to 
the Newtonian potential. (5$ is given by: 

The Fourier transform of 6{R,t)R is therefore: iaSk/k, 
and the Fourier transform of the DDE density contrast, 
6rf,, in the linear regime is given by: 



--{l + w)g{a) p 



where 



(19) 
(20) 



Motivated by observations, we assume (1 -I- w) to 
be small at late times, and (1 — ^m) to be small at 
early times. To a first approximation then we approx- 
imate /(a) by its ACDM value. At late times a very 
good approximation to f(a) in th e ACDM model is ^l'^ 
(|Peebles|[l980t iLahav et al.lll991[) . Note that corrections 
to this fitting formulae occ ur for quintessence- like models 
(jWang fc Steinhardtlll998f ). However, these are generally 
too small to greatly affect our results. Hence we take the 
ACDM expression for simplicity. Even when it is not ac- 
ceptable to approximate /(a) by its ACDM, we do not 
expect this approximation to greatly alter the qualitative 
nature of our results or the order of magnitude of 6^. 

Transforming back into real space we have: 



e(iin) (1 + w)g(a)»ma^g^ 



r'6{r',t)dr', (21) 



Our expression for (J^'^^ is independent of the mass of 
0. This is because we have assumed that ^ 0{H) 
in the linear regime and that the matter perturbation 
is small compared with H~^. The mass of the scalar 
field therefore has only a negligible effect on the scales 
over DDE clustering occurs. It is clear that the sign of 
6^k is the same as the sign of —g{a)dk- At late times 
g{a) > 0, and so a linear local overdensity of matter 
(0 < 5m ^ 1) corresponds to a DDE void. Similarly, 
there is a DDE overdensity at late times (in the linear 





|5j,/(Uw)| 


HR^(z=30)=0.03 ~~"^V^ 

8 (R=0, z=30) = 0.02 
m 
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Redshift: z 

X 10"^ Evolution of DDE density contrast in the linear regime 




Fig. 2. — Evolution of the DDE and matter density contrasts 
in the linear regime for an overdensity that begins to collapse at 
z = Zi = 30. We have taken Qrni^ = 0) = 0.27. Initially the 
matter density profile is taken to be Gaussian with Smizi,r = 
0) = 0.02 and HiVc = 0.03. This corresponds Smir = 0) = 0.3 and 
Rc = arc ~ 31/i~^Mpc today. \5^\ is always small compared to 
5m- At early times, the DDE density contrast is positive however it 
must be noted that our approximation is only valid in this case for 
z < 29. At late times it is negative. We can see that \S^/{1 + w)\ 
is growing at late times in the linear regime; indeed it is in fact 
growing faster than Sm for z < 0.23. 



regime) if there is a local void of matter. If an initial 
mean density perturbation dk{zi) begins to collapse at 
z — Zi when a — ai = 1, Qm = ^mi and Hi we have: 



1 



(22) 



317 



0.4 



0.2 



0.5 



a-^/']S,iz.), 



and so 



^<t>k 



9il + w)nm^Hf 



mi 



8 / ^7ni 

3 v r^7n 



4 

0.5 



Sk{zi) 



l + z 
1 + z, 



5/2N 



fc2 



Our expression for (50 was derived under the assumption 
that \Sij\R'^ < \Sv\R < \Sv\. As result it wiU break 
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down as we approach z = zi. More precisely, we find 
that our approximation is vahd at i? = for a > ai(l + 
2i?ji?_2(zj)) i.e. z < Zapp = Zi{l - 2HiR^2{zi)) where 
Hi is the initial value of 7? and R-2{zi) is the initial value 
of R-2- /.From this expression we can see that a matter 
overdensity corresponds initially to a DDE overdensity 

((J^''""* > 0), however this DDE overdensity becomes a 



DDE void at some z = Zcrit, and 6^ 
We see that Zcrit is given by: 



(lin) 



< for z < Zc 



1 + z, 



crit 



l + Zi 



(^crit) 



6 



3^l'^^{Zcrit)\ ^ 



where we have taken fir 



1. If QrniZcrit) = 1, which 



should be a good approximation for 1.8 < z < 3000, we 
have Zcrit — (3/8)^^^(1 + z^) — 1. For this analysis to be 
valid we must have Zcnt < -^app which with ^mizctit) ~ 1 
holds HiR-2{zi) ^ 0.24. For inhomogeneities that are 
larger than this at z = Zi we do not expect Zcrit to be an 
accurate approximation to the redshift when the DDE 
density contrast changes sign. 

Note that Zcrit does not depend on the size of the in- 
homogeneity, although, of course, we must have Zcrit 
Zapp for our analysis to be valid, which implies that at 
z = zcrit, HR_2 < 0.2. 
If ri„i = 1 exactly then (5^'" 



times. If 



< 1, then |(50/(1 
5c 



+ w) const at late 
+ w)\ grows likes: 



w 



cx 2r2; 



0.2 



0.2 



at late times in the linear regime. This implies that 
|(50/(1 + w)\ tends to a constant in the linear regime 
if rim — 1, and is growing for < 1. At late times 
in the hnear regime then: l^^/^ml oc Vl'^-^a^^. If we 

take r2„, = 0.27 today then IfJ^'""*/!! + is predicted 
to grow faster than 5m when z = 0.23. iDutta fc Maori 
(jlibO?) also found that at late times in the linear regime 
|(50/(l + w)| grows faster than 5m- This is behaviour were 
to continue into the non-linear regime it would, of course, 
imply that |(50/(1 + w)| would ultimately grow to be very 
large and indeed dominate over 5m- As we show in this 
paper, however, before this can happen the matter per- 
turbation goes non-linear and this slows the growth of 
\5l{l + w)\. 

We conclude our analysis of the linear regime by noting 
that Eq. |[2T|) implies that at late times: 



(23) 



(217jtf-an)H^(ar_2)'<5(r-2,i), 

where r_2 is defined the smallest value of r for which 
dhi5 / dlnr > —2, and Cs ~ 0{1) depends on the precise 
form of 6. 

We have shown that in the linear regime j^^j at the cen- 
tre of the inhomogeneity is always smaller than \6(r^2, t)\ 
by a factor of about (1 + w){Har^2)'^', Har^2 is roughly 
equal to the physical size of the inhomogeneity as a frac- 
tion of the horizon size. If the initial density perturbation 
is Gaussian i.e. 



5m{zi) = (5oexp(-r^/r; 



2\ 



then 



and: 



5.ir) 



^V^erf(^) - 2^e-(^/-)^) 



5T'\r = Q,t) 



9(1 + w)5QHfrl 
20 



(217, 



at late times. 

We plot the evolution of 5m{r ~ 0) and 5(j}{r = 0) 
for a Gaussian initial matter density perturbation in 
the linear regime in Figure [2] We have assumed that 
5m{r = 0,zO = 0.02 at Zi = 30, and H.r^ = 0.03, 
which corresponds to and 5m(r = 0, z = 0) = 0.3 and 
HaVc = 0.010 today i.e. Rc = ar^ ~ ilh'^ Mpc. For this 
choice of inhomogeneity, our approximati on for 5s is only 
oodfor 1-t-z < 0.94(H-Zi) w 29. Just as lDutta fc Maori 
200l saw in their numerical simulations, we find that 
close to z = Zi, (5^/(1 -I- w) > 0. It then becomes nega- 
tive, and at late times it continues to grow more negative, 
and it appears as if it might overtake 5m at some point 
in the future. However, as we shall show below, this does 
not occur and \5^\ remains small at all times. Note that 
the apparent increase in 5^ at early times is due in part 
to fact that our approximation breaks down as one ap- 
proaches z = Zi. If one wishes to study the evolution of 5^ 
all th e way back to z = Z j , the numerical approach taken 
bv iDutta fc Mao^ ([2007) provides accurate results. Our 
real focus in this work is, however, not on what occurs at 
very early times in the linear regime, but how the DDE 
density contrast evolves when the matter inhomogeneity 
goes non-linear. We consider this below. 



5.2. Dynamical Dark Energy Clustering in the 
Quasi-Linear Regime 

We begin our investigation of how 5^ evolves when the 
matter perturbation exits the linear regime by consid- 
ering a matter perturbation that is in tjie weakly non- 
linear or quasi-linear regime i.e. —\<5<1Q. Eqs.Q 
and (P give M{r,t)/R - M(r)/i? + 0((5|<5i, 5^51 525^). 
We suppose that at some initial time, ti, R t = HR and 
R = r. M(r) is then the mass inside the shell with radius 
r a,t t ~ tj. We write: 



M{r) 



Hy{l + 5,{r)). 



5i{r) is interpreted as the initial mean matter density 
contrast. In the full non-linear regime, which we con- 
sider in the next subsection, we are only able to evaluate 
R{r,t) and hence, via Eq. (flS)) . 5^, analytically in a 
matter-dominated background (r2,„ = 1). In many cases 
of interest, however, dm ~ e.g. superclusters and 

voids. Although the linear approximation fails for such 
objects, a good leading approximation to the true_mean 
density contrast of matter, 5, in the range, — 1 < (5 < 10 
is given by: 



l + 5{r,t) « (l-25ii„(r,i)/3) 



-3/2 



where is the linear mean density contrast. This 
approximation improves as flm and is exact for 
Qm = 0. When this approximation holds we say that 
we are in the quasi-linear regime. 
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Fig. 3. — Evolution of the DDE and matter density contrasts for 
a supercluster wliicli today has a core matter overdensity of 2 and 
a core radius of I5h~^ Mpc. We have take Q^n = 0.27 today. 15^ 
is always <SC 5^, and (5^/(1 + lo) < which corresponds to a DDE 
void. As in the linear regime, \5^/{l + grows monotonically 
with time, and at late times \5^/{l + w)\ is growing more quickly 
than 5m- 



Fig. 4. — Evolution of the DDE and matter density contrasts for 
a typical void which today has a core matter density contrast of 
—0.93 and a core radius of 15h~^ Mpc. We have taken Qrn = 0.27 
today. I (5^ I is always <^ \ Sm | , and 5^/ {l+w) > which corresponds 
to a DDE overdensity. As in the linear regime, 1 5^/(1 + w^)! grows 
monotonically with time and at late times \5^/{l + w)] is growing 
more quickly than 1 5m \ ■ 



In this section we find S^{S) for inhomogeneities in the 
quasi-Hnear regime. We are again assuming that any de- 
viations from the ACDM model for structure formation 
due the background DDE evolution are sub-leading or- 
der i.e. (1 — f2m)(l -I- w) is small. If this is not the case 
then our expression for in terms of d and Sv is still 
valid, but the evolution of S and Sv would change; even 
still the order of magnitude of S^j, would not be greatly 
effected. 

The peculiar velocity, 6v — R^t — HR, is related to 
S{r,t) thus: 



Sv 



ItR 



1 



HRf{am + 5f'^-l), 



where, as above, /(a) w at late times. In the quasi- 
linear regime the physical radius, R{r,t), of a shell with 
mass M{r) is given by: 



(24) 



Using Eq. (fT5|) , we find that the dark energy density con- 
trast, (5^, in the quasi- linear regime is well approximated 

by S^lf^ where: 



e(qi) 



(25) 



IR 



For 0.27 < < 1, the first term in Eq. (|25D is neg- 
ative definite for all positive values of_^ for which the 
quasi-linear approximation holds (0 < (5 < 10) and pos- 
itive definite for — 1 < ^ < 0. The second term in Eq. 
([25|) is clearly positive definite but vanishes at i? = 0. 
The relative magnitude of the two terms depends on the 
choice of initial density profile. Howejver, at i? = it is 
clear that an overdensity of matter, S > 0, corresponds 
to a dark energy void, S^{R — 0,t) < 0. Similarly, a 
void of matter corresponds to a local DDE overdensity 
at i? = 0. In Figures [3] and [4] respectively show the 
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Evolution of Perculiar Velocity in the Sptierical Collapse Model (SCM) 



10 



- Improved SCM 

- Standard SCM 




10 10 10 

Mean density contrast: 1+6 



Fig. 5. — Plot showing how the pecuHar velocity, hsc = 
1 — 5v/HR, varies with the mean density contrast, 5, in both the 
standard spherical collapse model (dashed line) and the improved 
one (solid line) of ! S haw Mota (2007). The improved model is a 
good fit to the results of N-body simulations for <5 > 15. We note 
that hsc reaches a maximum value of about 2 when 5 fii 18 in 
the improved model, whereas in the standard collapse model hgc 
continues to grow until S 140, where it reaches a maximum value 
of about 9.5. In the non-linear regime, the leading order expression 
for (5^ contains a term that is proportional to Sv^ oc h^^. There- 
fore if the standard spherical collapse model were used _to evaluate 
5^ it could result in it being over-estimated, for 20 < 5 < 140, by 
roughly 1-2 orders of magnitude. 



evolution of ^^/(l + w) and ^ at i? = for a typical su- 
percluster and a typical void. For both the supercluster 
and the void we assume an initial density contrast profile 
that is flat for r < Tc and drops as for r > Vc- We 
take the initial instant to be in the far past i.e. Zi — > oo. 
The supercluster is taken today to have a central density 
contrast of (5(i? = 0) = 2 and r = rc to corresponds to 
a physical radius of Rc = 15/i~^ Mpc, which are appro- 
priate for an object such as the local super-cluster. For 
the void we take typical values of = 0) = —0.93 
today and Rc — 15/i~^Mpc. In each case we have 
taken ^rn{z = 0) = 0.27. It is clear that in both cases 
|<5</)/(l + w)\ grows faster than \5m\^ however the DDE 
density contrast remains much smaller than \5m\ at all 
times. Note that 5^ evolves monotonically for objects in 
both the linear and quasi- linear regimes. 

5.3. Dynamical Dark Energy Clustering in the 
Non-Linear Regime 

The quasi-linear approximation to 5^ is valid for — 1 
5 < 10, however it breaks down for large positive values 
of the matter density contrast, 5. To go further we evalu- 
ate 5^ in the fully non-linear regime. To do this we need 
to calculate both the evolution of the peculiar velocity, 
&v, and that of the mean density contrast mass contrast, 
6. The non-linear evolution of the matter overdensity 
is significantly more simple to calculate analytically in 
a matter dominated Universe {p.m — 1)- Fortunately, 
when the matter density contrast is large today 5 > 100, 
the evolution of both R t and e^™' in the inhomogeneous 
region are well-approximated by the f2„j = 1 solution. 
This is because the epoch of dark energy domination be- 
gins sometime after the cluster has turned around. 



If VLr a = 1, we are effectively d ealing with a Tolman- 
Bondi (|Tolmanlll93l lBondilll947l ) background at leading 
order. We therefore have analytical solutions for R and 
t in terms of r and a parameter 77 (r, t): 

R^^M^r [l-cos^), (26) 
25i{r) 

t = ts{r) + ^^^§^{ri-^inr,), (27) 



where 



2H4'\r) 



o{st) 



and 



l-5,ir) 
l + S,ir) 



25]'\r)~-jf\r)+0{5l/\ 



We assume that at early times, (z > 1.8), « 1. These 
solutions then provide an excellent approximation to the 
true early time evolution of the inhomogeneity. These so- 
lutions also provide an good approximation at late times, 
provided we are dealing with an overdensity of matter 
that turns around prior to the onset of dark energy dom- 
ination {5 > 100 today). Assuming that we are dealing 
with such an overdensity, we find: 



5M 



2 



H^R 



9ia/aif{r] - sinry)^ 
2(iA»)2(l-<5t)2(l- cos 77)3 



1 



1 



(29) 



IS 



1. 



matter-dominated , 
In the non-linear 



When the background 
{a/ai)^ /{t/t,f = 1 and Q.„^ 
regime 5t = ts{r)/t <^ 1, and terms proportional to 6t in 
the above expression can be dropped. 

We have so far modelled the matter as being a spheri- 
cally symmetric, pressureless perfect fluid. This approx- 
imation neglects the random motion of the matter par- 
ticles, and other effects due to the break down of spher- 
ical symmetry. These are small prior to turnaround but 
thereafter act to slow down the collapse of the inhomo- 
geneity. Without any such corrections the matter pertur- 
bation would eventually collapse to a point. When these 
effects are included however the inhomogeneity virialises 
and relaxes to a steady-state. In a matter-dominated 
background, each shell of constant M ceases to collapse 
and virialises when its radius, i?vir, is about one half of 
the radius at which it turned around, i?ta- N-bod y sim u- 
lation s suggest that i?ta/-Rvir ~ 1-8 (see (Hamilton et all 
[1991 . 

Since process of virialisation is not accounted for in our 
matter model, we must account for it in some fashion 
if we are to make accurate predictions. This is gener- 
ally done via the rather ad hoc process of manually halt- 
ing the collapse when each shell reaches its virial radius, 
i?vir; in this model one takes i?vir = Rta/2; this occur at 
77 = 37r/2. Whilst this provides a reasonably good ap- 
proximation for the magnitude of the density contrast, 
S, the peculiar velocity, Sv, is discontinuous at instant of 
virialisation (77 — 37r/2). This ad hoc procedure there- 
fore gives widely inaccurate predictions for the peculiar 
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Growth of Density Contrast with the Linear Density Contrast 




2 3 4 

Linear Density Contrast: 5 



Fig. 6. — Plot showing how the true mean density contrast, S, 
depends on the hnear one, <5iin oc a{t) for Qm = 1- We show both 
the standard spherical collapse m odel (SCM), th e dashed line, and 
the improved SCM proposed by IShaw fc Motal (2007), the solid 
line. For 15 < 5 < 140, the S grows faster in the standard SCM 
than in the improved one, however the overall discrepancy between 
the models' predictions for 5((5iin) is not nearly as great as it is for 
their predictions for hsci^) = 1 ~ Sv/HR: see FIG. [5] 



velocity, Sv = Rt — HR (see Figure[S]). Since our formula 
for depends on both Sv and S we require a more realis- 
tic model for virialisation in the context of spherical col- 
lapse. Such an impr oved model was proposed recently by 
IShaw fc Motal (1200711 . In this improved spherical collapse 
model, it is assumed that hsc = Sv/HR is a function of 
S i.e. hsc = hsciS)- R is still parametrized according to 
Eq. l[26|l . and it is found that: 



2H.,S^'\r) 



Tir), 



where t = rj — sin 77. The function T{t) is then fixed so 
that: Prior to turnaround, when deviations from spher- 
ical symmetry are small, T(t) r. For larger values 
of S (> 15), T{t) is chosen so that hsc{S) agrees, un- 
der certain reasonable and common assumptions, with 
the behaviour that is seen in N-body sim ulations of 
iHamilton et~aI1 (|1991f ) . IShaw fc Motal (|2007t l provide the 
following fitting formula for T{r]): 



T{ti{t)) = t + 



3.468(r/ - t)-i/2 exp (-ifcll) 



(1 + 0.8(r/ - t)i/2 _ 0.4(r/ - r)) ' 

(30) 

where r/ = 5.516. In this model Rvh/Rta. ~ 0.59 
which compares more favourably wi th the value of 0.56 
suggested by N-body simulations ( Hamilton et al.|[l991f ) 
than the value of 0.50 that is generally used. In this im- 
proved model, S^ is still given by Eq. [HI and for n„i = 1 
the peculiar velocity is given by: 



St 



R.t 



- 1 = 



1 



sm?7 



dr 



HR HR ' V 2 (l-cosr;)i/2dr 
and the density contrast is given by: 

9r(r/)2 



1, 



1 + S = 



Figure [5| shows hsc = 1 — Sv/HR as a function of S 
in both the improved spherical collapse model, and the 
standard one (with virialisation occurring suddenly at 
R = i?vir)- Note that in the standard spherical collapse 
model (SCM), hsc continues to grow until 5 « 140, and 
then drops sharply, whereas_in the improved model hsc 
reaches a maximum when (5 « 18. Note also that the 
maximum value of ft.sc predicted by standard SCM is 
almost five times as large as the maximum value gives 
by the improved model. Since S^ depends in part on Sv"^ 
and hence h'^Q, using the standard SCM instead over the 
improved model could lead to S^ being over estimated by 
as much as 2500%. 

Figure [6| shows 5, for — 1, as a function of the mean 
linear density contrast, Sw^ oc a, in both the improved 
model and the standard one with virialisation put in by 
hand. We use our modified spherical collapse model to 
find the evolution of (5^. We find then that the DDE 
density contrast in the non-linear regime is given by: 

^k£!ll±i^ r -S[R' ,t)G{5-n^)R' dR' 



+H\1 



where 



GiS;n^) 



w)R^h{Sf, 
6h{S) 



(31) 



1. 



5 

When (5 < 1, G{S) 1 - 2{n~^°-\- 1), and deep in the 

non- linear regime, when (5 3> 1, G{S) ^ 1. In the flm — 1 
case we find that the following fitting formula is accurate 
to within 2%: 

G{S; 17„ = 1, a = 1) w tanh (-0.9218-h (32) 

0.543^6 \ _ 0.38 sin^ (ln(l) - 2.49) \ 
403 + J (ln((5) - 2.49)3 j ' 



195 



2(1 — cos 77)3 



We note that when — 1, G{S) < for S < 10.6 and is 
otherwise positive. In a cluster with a core overdensity 
> 20, 5^ is therefore generally positive at i? = i.e. there 
is DDE overdensity at i? = 0, whereas for clusters with a 
core overdensity < 10 — 20, there is generally a DDE void 
around R = 0. Even in the deep non-linear regime, the 
central DDE overdensity is still surrounded by a void, 
with 5^ becoming negative when 5{R) « 10 — 20. We 
see that \S^\ - 0{n^{l + w)H'^Rl5/2), where Rc is the 
radius of the core of cluster (roughly defined by smallest 
radius after which S decreases faster than 1/R^), and 5c 
is the core density contrast i.e. the density contrast at 
R = Rc. 

We use Eq. ([5T|l to plot the evolution of (50 (i? = 0) in 
the non-linear regime. For 5^1, the matter overdensity 
in a background with 0.2 < Vim < 1 evolves, to a good 
approximation, according to: 

l + S{7j){n,n)^F{a){l + SQ{v)), 
where Sq = 9T'^/2{1 - cosrjf and F{a) = g^i^^^^H-^ ~ 
f7-0-4^ We expect the largest contributions to S^{R = 0) 
to come from regions where 5 3> 1 which implies that the 
we expect G{S; ilm) ~ 1 for all flm- Therefore, using Eq. 
(j3ip . we approximate: 



So{v)G{So{v),^m = l)R'ir,v)dR'. 
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Fig. 7. — Evolution of the DDE and matter density contrasts for 
a cluster which today has a core matter overdensity of 200 and a 
core radius of 1.5/i~ Mpc. We have taken = 0.27 today. 15^ 
is always <C Sm- At late times S^/(l + w) > 0, which corresponds 
to a central DDE overdensity, however at early times it is negative 
and there is a DDE void. We see that changes sign when the 
mean core overdensity of matter is ~ 0(10). 



Fig. 8. — Evolution of the DDE and matter density contrasts for 
a cluster which today has a core matter overdensity of 1.6 X 10'* 
and a core radius of 0.18/i~^ Mpc. We have taken = 0.27 
today. \S^\ is always <^ 5m- At late times S^/{1 + w) > 0, which 
corresponds to a central DDE overdensity, however at early times 
it is negative and there is a DDE void. We see that changes 
sign when the mean core overdensity of matter is ~ 0(10). 



For simplicity, we take the initial density profile to have 
a central core with Si{r) — Si for r < Vc and for r > Tc, 
we take: Si oc r~^. 

In Figure[7]we plot S{R = 0) and ^^(i? = 0) for a clus- 
ter where today r = corresponds to a physical radius of 
1.5h~^Mpc, and the mean density contrast in R < Rc to- 
day is 200. In Figure. [8]we plot d{R = 0) and 5^{R = 0) 
for a_smaller but denser cluster with Rc = 0.18/i~^Mpc 
and 5{R = 0) = 1.6 x 10^ today. These values are fairly 
typical of clusters such as Coma which has been esti- 
mated to have a mean density contrast of 200 inside a 
radius of 1.5/i^^Mpc, and a core mean density contrast 
of about 1.6 X 10^ inside a radius of 0.18/i~^ Mpc. We can 
see that in both cases ^^/(l + w)^ 0(10^'"') today and 
that, as expected, it is positive. In both cases, 5^/{l + w) 
was negative in the past, changing sign when the mean 
overdensity of the core was 5 ^ 0(10). The core in these 
cases corresponds to the flat part of 5 i.e. R < Rc- For 
more general density profiles, the core radius would be 
taken to correspond to i? < i?_2 where i?_2 is the small- 
est value of R for which 6 faster as R^^. In all of the 
plots we have taken ilm{z = 0) = 0.27. 



We can see that in both cases S^/{1 + w) continues 
to increase today. This increase is due to the continued 
collapse of matter onto the collapsed core. At very late 
times, however, this accretion will cease, and i5^(i? — 
0)/{l + w) const for all ri,„. We can see this from Eq. 
(j3ip . At late times, and assuming minimal accretion: 
S — a'^A{R), for some function A(i?), and so: 

„3q /72 poo 

S^/a + t / R'W dR', (34) 

and a^flmH^ = const which implies 15^/(1 + w) ^ const 
at late times. It is clear that S^/{1 + w) is smaller than 
(5-2 — 5{R — R-2,t) by at a factor of about QmH'^R-2- 
In follows that even at late times 6^/{l + w) <C 1, 
as 5^/{l + w) ^ 0{nmH^Rl^S^2) 0(GM_2/i?-2), 
where A/_2 ^ M{R^ i?_2). 

5.4. Discussion 

We have seen that in all regimes (linear, quasi-linear 
and fully non-linear regimes), the magnitude of DDE 
density contrast scales as: ilmH^Rc^c, where Rc is the 
radial scale of the cluster, and Sc is the mean matter 
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density contrast m R < Re- Rc is roughly defined to 
be equal to R^2, which is the smahest radius for which 
5c first decreases faster than 1/R^. If an inhomogencity 
is sub-horizon when it begins to collapse, we found that 
shorter after the collapse begins 6 and have the same 
sign i.e. a matter overdensity equates to a DDE over- 
density, and the same for voids. However, a certain time 
after collapse begins, changes sign. This change of 
sign is related to the fact that at late times the decaying 
mode of the matter perturbation is negligible. When this 
occurs, a mean over density of matter results in a DDE 
void, and vice versa. Deep in the non-linear regime, when 
(5c > 10 — 20, however, 6^{R — 0) changes sign again. A 
large mean overdensity of matter then corresponds to an 
overdensity of dark energy at i? = 0. 
In both the linear and the quasi-linear regimes, | (50/(1-1- 
increases faster than 6c, although it is always much 
sm aller than Sc- This behaviour was also observed 
by iDutta fc Maori (|2007[ ) in their numerical study of 
the linear regime, and it lead them to wonder whether 
\6^/{l + w)\ might grow to be 0(1) or larger in the non- 
linear regime. Using our results, we have shown that the 
onset of non-linear evolution for the matter perturbation 
slows down the growth of \d^/{l + w)\ and that at very 
late times S^/{1 + w) ^ const and remains <C 1. In the 
next section we consider the profile of S^/{1 + w) in and 
around typical astrophysical objects. 

6. ASTROPHYSICAL APPLICATIONS 

In this section we use our results to evaluate the profile 
of the DDE density contrast around typical clusters, su- 
perclusters and voids of matter. We also investigate the 
possibility of detecting the dark energy clustering using 
the Integrated Sachs Wolfe (ISW) effect. 

6.1. Galaxy Clusters 

We begin by considering the DDE density contrast 
in a virialised gal axy cluster with similar proper ties to 
the Coma cluster. iNavarro. Frenk fc White! (|1997D (here- 
after NEW) used high resolution N-body simulations to 
derive a universal density profile for the dark-matter ha- 
los of virialised structures such as galaxies and clusters. 
The resulting NEW profile is given by: 



iR/Rc)il + iR/RcW 

where Rc defines the scale at which 6m changes from a 
1/r drop-off to a one, and the characteristic density 
contrast 6c is given by: 

_ 200 c3 

' ^ ~3~ln(l + c) -c/(l + c)' 

where c — R200/ Rc defines the concentration of the halo, 
and i?2oo is roughly the virial radius and is defined to be 
the largest value of R for which the mean density contrast 
is > 200. The mean density profile, (5, is then given by: 

- _ 3Je(ln(l + R/Rc) ~ Rc/{1 + R/Rc)) 
{R/Rcf 

Even though 5 cx 1/i? for R < Rc, the integrands in Eq. 
([5T|) still behave well enough for small R for our analysis 
and our expression for 6^ to be valid and applicable. The 
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Fig. 9. — The DDE (a) and matter density (b) profiles for a 
Coma-like cluster with NFW profile: Rc = 0.18/i~^Mpc, R200 = 
1.5/i~^Mpc. We have taken flrn = 0.27. See text for discussion. 



DDE density contrast in a virialised cluster is given by 
Eq. ([3T|) . We note that since 6v'^ (x h (x flll"^, in the 
non-linear regime, and 6 depends on R only though R/Rc 
t hat, for fixed 6c, 6^ (x f^m/ l + w)H'^Rc- 

iGeller. Diaferio fc Kurtg £1999) studied the mass pro- 
file of the Coma galaxy cluster and found that over the 
entire range of i? < 10Mpc/i~^ it increased with r at 
the rate predicted by the NFW profile. They found that 
three different observational samples were well fitted by 
Rc = 0.182 ± 0.030, Rc = 0.167 ± 0.029 and 0.192 ± 
0.035/i~^Mpc with, in all cases, -R200 ~ 1.5/i^^Mpc. In 
what follows we take Rc to be the mean of three values, 
Rc « O.lSft."^ Mpc, and i?200 = 1.5/1"^ Mpc, which cor- 
responds to c « 8.3. We plot the profile of DDE density 
contrast for a Coma-like cluster in Figure [51 R = Rc and 
R = ^200 arc indicated on the plo ts. We have taken 
= 0.27 in line with WMAP (Spergel et all l2003l . 
|2007| ). We note that 6^ is positive around R = and 
also that it is very small 6rj,{R = 0) « 3 x 10~^. As R 
increases, (5^ becomes negative when R « 1.5i?200) and 
for R > 4i?2oo it tends to from below as i? — > 00 in line 
with our expectations. 

6.2. Superclusters 
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Superclusters are the largest known gravitationally 
bound massive structures. The dark matter halos of su- 
perclusters typically have radii of about 10 — 25h~^ Mpc 
and density contrasts of about 1 ^ (5 < 15. We con- 
sider two examples of such objects: the local superclus- 
ter (LSC), of which our galaxy is a part, and the Shapley 
supercluster (SSC) which lies about 650 Mlys away. The 
LSC has a mean overdensity of 5 y 2 — 3 over a scale 
- 15/1-1 Mpc (Hoffman 1989: Iruilvl [198^ . The SSC has 
been found to hav e an overdensity of S ^ 1 0.3 over a scale 
of lO lfe-^Mpc (iBardein et all l2000( l. iBardelli et all 
(|2000( ) also found that if the SSC had evolved linear then 
it would today have a linear overdensity of Sssciin ~ 1-3 
and a linear scale of « 21.6/1-^ Mpc. For simplicity we 
model the matter density contrast of the clusters thus: 

• The clusters have a homogeneous core with physi- 
cal radius Rc- In R < Rc, S = Sq and Sun = ^Uno- 

• For R > Rc, the initial density contrast drops off 
as where R = Rcx{l + 5o)/(l + S{R)Y'^. 

With this choice of density contrast, the evolution of the 
cluster is free from both shell-crossing and shell-focusing 
singularities in the_ linear and quasi-linear regimes. For 
the LSC we take 5q = 2 and Rc = lbh~^ Mpc and for 
the SSC we take 5o = 10 and Rc = lO/i"! Mpc. We 
additionally considered the S300 st ructure, in the vicinit y 
of the SSC, which was identified bv lBardelli et all (|2000D . 
This structure has a mean overdensity o f 5 1.9 on a 
scale of 24.8/1-1 Mpc (iBardelh et al.ll200G[ ): we take 5q = 
2, Rc = 25h~^ Mpc for this structure. 

We plot the DDE density contrast profile for the LSC, 
SSC and S300 structures in Figure [TO] and Figure [TTl 
The shape of 6^ is similar for all three objects, and 
the magnitude of the DDE density contrast scales as 
(1 + w)5qH^ R^\ it is largest for the S300 structure and 
smallest for the LSC. In all three cases, it is clear that 
1(5^/(1 + w)] is small (< 6.2 x 10^^) and negative which 
corresponds to a local DDE void, 

6.3. Voids of Matter 

In addition to local overdensities of matter, the Uni- 
verse also contains localized underdensities of matter or 
voids. Voids typically have the similar radii to super- 
clusters, although they can also be much larger. Since 
1(50/(1 + w)\ is proportional to the square of the radius 
of the inhomogeneity, we expect voids to induce some 
of the largest DDE inhomogeneities. The classic histor- 
ical example of a voi d is also one of the larg est and was 
discovered in 1981 bv lKirschner et al] (|1981l ) in the con- 
stellations of Bootes and Corona Borealis. Recent mea- 
surements of this Bootes void hav e found it to be roughl y 
spherical with radius 62h^^ Mpc (jKirschner et al.l[l987f ). 
21 galaxies have been observed in the Bootes void and 
ij;s mean de nsity contrast is estimated to be: —0.80 < 
6 < -0.66 (Dev. S trauss fc Huchral[l990[ ). The Bootes 
void is particularly large and is sometimes termed a su- 
pervoid. Data from the 2dFGRS shows that the aver- 
age radii of voids in NGP (North Galactic Pole) and in 
SGP (South Galactic Pole) are 1 4.89 ± 2.67/t-i Mpc and 
15.61±2.48/i-iMpc respectivelv jHovle fc Vogelevll200 j: 
iBoleiko. Krasihski fc Hellabv .2001 . The average mean 
density contrast for these voids is 6 = —0.94 ± 0.02 in 
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Fig. 10. — The profiles of DDE and matter density contrasts for 
superclusters. We have taken fl^ = 0.27. Figures (a) the DDE for 
the local supercluster (LSC) and Figure (b) shows our model for 
the LSC's matter density profile. We have taken the LSC to have 
a mean matter density contrast of 2 inside a radius of I5h~^ Mpc. 



NGP and S = -0.93 ± 0.02 in SGP (|Hovle fc VogelevI 
[200l . 

Despite the large radii of these voids and supervoids, 
their scales are still very much sub-horizon and, as such, 
our results are applicable to them. The evolution of 
voids is well approximated by the quasi-linear regime 
and so we evaluate using Eq. We model 

both voids and supervoids as_ having a Gaussian ini- 
tial mean density profile i.e. Sun oc exp(— a;^/2) where 
R^ Rcx{1 + Sq)/{1+5{R))^/^ defines a;. With this choice 
the evolution of the void is free frorn both shell-crossing 
and focusing singularities provided 6o > —0.95. 

In Figure [1^ we show S^j, and 5 for both a super- 
void with similar properties to the Bootes void {Rc = 
62/i-iMpc, So = -0.75), and Figure [13] shows the 
same thing but for an average void {Rc — 15/1^^ Mpc, 
Sq = —0.93). In both cases we see, as expected, that 
(5^ > and that a void in the matter distribution cor- 
responds to a DDE overdensity. We note that the mag- 
nitude of \S^/ {1 + w)\ is an order of magnitude larger 
for a 60h~^ Mpc radius supervoid than for an average 
15/i-iMpc void. At the centre of a supervoid we see 
that (5^ ~ 9 X 10-^(1 + w), whereas at the centre of an 
average void (5^ w 5 x 10^^(1 + w). In all cases is both 
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Fig. 11. — The profiles of DDE and matter density contrasts 
for the Shapley supercluster (SSC). We have taken = 0.27. 
Figure (a) shows the DDE profile for the Shapley supercluster, 
which has been modelled as having a mean density contrast of 10 
inside a radius of 10/i~^ Mpc. Figure (b) shows the DDE profile 
the S300 structure (with mean density contrast of 2 inside a radius 
of 25/1-1 Mpc). 



positive and very small. 

Recently, there has been a great deal of discussion, 
both theoretical and observational, about the possibility 
that extremely large voids ~ 100— SOO/t"^ Mp c ) exist 
within the visib le Universe (Rudnick. Brown fc Williamd 
I2007t llnoue fc Silk 2006 ). Additionally, a number of au- 
thors have spe culated that we may even be living inside 
such an ob ject ("Schwarz fc Wcinhorst"2007':'C onlev et all 
I2007t iZeha vi ct al. 1998: GiovancUi ct al. 1999|). For ex- 
ample, ^Iidnick, Brown fc Williams (2007) showed that 
both the WMAP cold spot and an observed dip in the 
surface brightness and number counts of extragalactic 
radio sources could be explained by the presence of an 
almost comp l etely empty Rc « 105/1"^ Mpc at z < 1. 
llnoue fc SiiM (|2006l ) showed that the ISW effect due to 
a void with (5o = -0.3 and = 200 - imhr^ Mpc 
would be observed as cold spots with a temperature 
anisotropy AT/T ~ 0(10"^), which might explain the 
observed large-angle CMB anomalies. If extremely large 
voids such as these exist they would be associated with 
the largest DDE overdensities within the visible Uni- 
verse. In Figure [14] we plot the 5^ profile for both 



Fig. 12. — The DDE (a) and matter (b) density contrast profiles 
for a supervoid. We taken Cm = 0.27 and modelled the supervoid 
on the Bootes void: <5(0) = -0.75 and Rc = 62/1"! Mpc. See text 
for further discussion. 



the Rc = 105/i ^Mpc, i5n ~ — 1 o bject postulated by 
iRudnick. Brown fc Williamsl (|2007t ) and f or the Rc = 
300fe- ^ Mpc, (5o ~ —0.3 void considered bv llnoue fc Silkl 
(|2006[ ). In the former case we find (5^(7? = 0)/(l + w) « 

2.4 X 10"^ and in the latter S^(R = 0)/(l + w) « 

9.5 X 10^^. These values of S^/{1 + w) are about two 
orders of magnitude larger than those associated with a 
typical void. However, even around such extremely large 
voids we still find that the DDE perturbation is small: 

+ ^ 10"'- 

6.4. ISW Effect due to DDE clustering 

We have seen that the DDE density contrast, 5^, is 
always small with magnitude < lO^'^ for average sized 
structures in the Universe. The Integrated Sachs Wolfe 
(ISW) effect is highly sensitive to the presence of dark 
energy. In a purely matter dominated universe, the ISW 
effect vanishes entirely for linear perturbations. We shall 
see that the DDE contribution to the ISW from a cluster 
with radius Rc is proportional to {HRc)^ and so we con- 
sider only large (but still sub-horizon) scales. On such 
scales the matter perturbation is small and we are in the 
linear regime. The ISW effect is most easily derived in 
the conformal Newtonian gauge in which the metric takes 
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Fig. 13. — The DDE (a) and matter (b) density contrast profiles 
for a typical void. \Ve have taken Qm = 0.27 and the void to have 
typical properties: <5(0) = —0.93 and Rc = 15h~^ Mpc. See text 
for further discussion. 



Fig. 14. — The DDE density contrast for extremely large voids. 
In FIG. (a) we consider a 105h~^ Mpc void with (5(0) = -0.95 and 
in FIG. (b) we have taken void to have radius 300/i~^ Mpc and 
central density contrast 5(0) = —0.3. The existence of voids with 
these properties as recently been postulated as a mean to explain 
the observed large-angle CMB anomalies. 



the form: 

ds^ = [(1 + 2*) dr^ - (1 - 2$) dx^] . 

Since we are dealing with shearless matter, to leading 
order, \1/ = $ and, in momentum space: 

2 Sftm'H^ r 3(1 — flrn)'H'^ 

k $ = d,„ 



5 a + 



2 2 
3(l + w;)H2/(a)'5„ 



(35) 



where Ti = aH — dlna/dr and 6m and 5^ are defined 
in the rest frame of the matter particles. Using Eq. (fTT" 
we see that: 



,2^ 3n„,n\ 9{i + w){i ^ n^jn^H'^Sr, 

k 9 = d„, 



2 



4fc2 



(36) 



The temperature anisotropy due to the ISW effect is 
given by: 

/■tq 

Siswin) = 2 $'((to - T)n, t) dr, 

Tt. Q 



where tq is the value of r today and tls is its value 
on the surface of last scattering and $' — d^/drj. The 
contribution to $' from the matter perturbation is: 

C = ^ (l-/(a))^„„ 

and /(a) « fi^^ at late times, whereas the contribution 
to <i>' from the DDE clustering is: 

- - ^MnUHm ^ ^^(2 + 3z« - /(a)) - w^)) 



Ak^ 

inmH\i - f{a)) 

2fc2 



v{n 



where Wp — dw / d In a and 



V{Qm,W,Wp) 



(1 - a„)(2 + 3w- f{a) - wp/il + w)) 



(2/(a)-r!™)(l-/(a)) 
It follows that: 

% 3(1 + w)H'^ 



2fc2 

(l + f^r + Y^-3(l + z.)), 
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assuming that Wp/{l + w) ^ 0{1) or smaller in the recent 
past then: 



< 



We found above that the largest DDE perturbations 
would mostly likely be associated with large voids of mat- 
ter, with « 9.5 X 10-*(1 + w) for a 300 Mpc void with 
central matte r density contrast = — 0.3 (such as that 
considered by llnoue" k SiM20M\ . In this case we would 
have: 



« -7 X lO^-'fl + w) 1.5 - 3(1 + w) + ^ 

*m V 1 + 

which is at most ©(— 10~^(1 + w)). In the absence of 
DDE clustering, the ISW temperature anisotropy of such 
a void was found to be \ AT/T\ ~ 0(10-^). We therefore 
predict that the DDE clustering contribution to the tem- 
perature anisotropy would be at most 0(10^^(1 -I- w)). 
For smaller objects (with radius Rc) the DDE ISW tem- 
perature anisotropy scales as R'^. 

In conclusion: the ISW due to sub-horizon DDE clus- 
tering in models where the dark energy is due to a mini- 
mally coupled scalar field with canonical kinetic term is 
always suppressed by a factor of about (1 -t- w){HRc)^ 
relative to the ISW effect to the matter perturbation it- 
self. It therefore seems unlikely that the DDE contribu- 
tion detected against the background of the ISW effect 
caused by the matter. 

6.5. Discussion 

In this section we have applied our results to evalu- 
ate the profile and evolution of ^^/(l -I- w) in the vicin- 
ity of typical galaxy clusters, superclusters and voids. 
We found that, today, the largest inhomogeneities in the 
DDE energy density value should occu r in the vicinity 
of sup ervoids such as the one in Bootes (jKirschner et al.l 
Il981| ). At the centre of a supervoid we predict 5^/(1 -f 
w) K, 1.6 X 10"''. At the centre of a galaxy cluster, su- 
percluster or average sized void, we found |(50/(1 -I- w)| « 
1 — 6 X 10^^; with (5^ > for both virialised clusters 
and voids, and 5^ < Q for superclusters. Generally 
\5^/{l + w)\ < 10^^ for supervoids, and < 1 — 6x 10^^ for 
other structures. Recently there has been a fair amount 
of interest in the possibility that extremely large voids 
might exists with radii Rc = 100 — 300 Mpc. Even 
for objects as large as this, however, we found that 
S^/{l + w) < 10-3. 

Given the current bounds on the dark energy (EoS) 
parameter: w{z < 1) = —1 ± 0.1, our results imply that 
if dark energy is described by a canonical scalar field that 
does not couple to (or only couples very weakly to) mat- 
ter and is minimally coupled to gravity, then its energy 
density today is virtually homogeneous with largest fluc- 
tuations no larger than 1 part in 10^ (or possibly one 
part in 10"* if 300 Mpc voids exist). Indeed, typical fluc- 
tuations in the DDE energy density will be smaller than 
a few parts in 10^. We found that superclusters are gen- 
erally associated with small, local DDE voids whereas 
most other structures (galaxies and voids of matter) are 
associated with local overdensities of dark energy. 

7. SUMMARY AND CONCLUSIONS 



In this paper we have studied how dynamical dark en- 
ergy clusters over sub-horizon scales in theories where the 
dark energy is described by a canonical scalar field with 
no coupling to matter. Sub-horizon dark energy clus- 
tering in such models has always been expected to be 
small, but previous studies of it have assumed the mat- 
ter perturbation to be small and have, for the most part 
relied on numerical simulations. In this paper, however, 
we used the powerful method of matched asymptotic ex- 
pansions to derive an analytic expression, Eq. (|18p . for 
the DDE density contrast in terms of the peculiar veloc- 
ity and density contrast of the matter. Our equation for 
holds not only in the linear regime, where the mat- 
ter density contrast is small, (5m ^ 1, but also in the 
far more interesting but far less studied quasi-linear and 
non- linear regimes where Sm ~ 0{1) or greater. 

There were then essentially two parts to this paper, 
finding an expression for 6^ as a function of Sv^ 6m, and 
H, and then evaluating this. The evaluating part is only 
done approximately for models that look sufficiently like 
ACDM. So that to a first approximation we can use 
ACDM to get H, Sv and Sm- Different dark energy 
models would give different predictions but just from the 
form of Scj) one can see that unless very large deviations 
from the ACDM predictions of Sv, 6m and H occur then 
Scji will be well approximated by taking a ACDM back- 
ground. 

Our key assumption in deriving our formula for 5^ was 
that the matter density contrast is only C(l) or greater 
on scales that are small compared to the Hubble length, 
H^^. This assumption holds very well for typical struc- 
tures such as galaxies, clusters, superclusters and voids. 
We also made a number of simplifying assumptions: 

• We assumed spherical symmetry. Although our 
interest is in dark energy rather than varying- 
constants, the sub-horizon dynamics of perturba- 
tions in the DDE scalar field are very similar to 
those of the scalar fields in varying-constant mod- 
els provided one sets the matter couplin g in such 
theories to zero. (|Shaw fc Barrowll2006allbi B stud- 
ied the evolution of a light scalar field in an inho- 
mogeneous background in the context of varying- 
constant. Both spherically symmetric backgrounds 
and ones without any symmetries were considered, 
and it was found that deviations from spherically 
symmetry did not alter the leading order behaviour 
of the scalar field perturbations when the matter 
coupling vanishes. Therefore, whilst the assump- 
tion of spherical symmetry can be relaxed, we do 
not expect it to greatly alter our results. 

• We assumed that we were not near a black hole 
horizon. This assumption allowed us to take 
CM/R ^ 1 on sub-horizon scales. This assump- 
tion could also be relaxed, and once again both the 
form of 6cf, and its magnitude are not expected to 
be greatly effected. We again base this expectation 
on the results found for light scalars in the con text 
of varying-constants bv I Barrow" k Shawl ((200l . 

• We assumed that the mass of the scalar field, m, 
was very small compared to l/Rc, where Rc in the 
length scale of the peak in the matter density con- 
trast. This assumption could be done away with. 
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however it would significantly complicate the anal- 
ysis in the non-linear regime. Moreover since we 
have require HRc <^ 1 and in most quintessence 
models m ^ 0{H), it is for most models gener- 
ally true that mRc ^ 1. If mRc ^ C(l) or larger, 
then we expect that it will result in a smaller DDE 
density contrast, i.e. smaller values of \S^\. 

One of the main motivatio n s for t his work was the re- 
cent work of iDutta fc Maori (|2007t l. They studied the 
clustering of a subclass of the DDE models which we 
have looked at in this work. They numerically evolved 
the linearized field equations for both (jj and the metric. 
In the regime where the linearized equations are valid 
(i.e. Sm <SC 1), our analytical results agree with theirs: 
shortly after the collapse begins, 6^/{l + w) is positive, 
but then as the dying mode of an overdensity of matter 
becomes small, S^/{l+w) changes sign and becomes neg- 
ative. At late times in the linear regime an overdensity 
of matter corresponds to a DDE void, and vice versa. 
We saw in this paper that this correspondence contin- 
ues to hold in the quasi-linear regime, — 1 ^ <5jn ^ 10, 
that is appropriate for superclusters and voids of matter. 
When dm ~ C'(l), Dutta and Maor found that \S^\ could 
grow to be relatively large (of the order of 0.01 — 0.1), 
and that this could lead to corrections to the equation of 
state parameter, of the order of a few percent. How- 
ever, as we noted in the introduction, the C(l) values of 
S„i for which this effect was found fall outside the regime 
where the linearized field equations which they used are 
actually valid. One should point out however, that in 
a matter dominated universe (flm ~ 1)? there is a clear 
mapping between the the linear regime of an overdensity 
and the exact fully developed non- linear overdensity. 
Presumably, such mapping should occur also in the pres- 
ence of a scalar field, though the multiple change of sign 
of the DDE perturbation complicate things considerably 
and this relation is not straight forward. We also saw 
that at late times in the linear and quasi-linear regimes 
I (50/(1 -I- grows faster than \6m\- If it were not for 
the onset of truly non-linear behaviour in the evolution 
of S,n, this woul d indeed lead to large values of l^^l seen 
at late times bv lDutta fc Maoil (|2007D . 

In our study we found that 5^/(1 -I- w) for reahstic 
astrophysical objects is very small, indee d it is generally 
<OnO-'^). The fast growth in seen bv lDutta fc Maoj 
(|2007) for Srn ^ 0(1) should therefore been interpreted 
as resulting from applying linearized field equations out- 
side their realm of validity, rather than a true physical 
effect. Furthermore, at very late times when the matter 
perturbation has virialised and accretion onto the core 
has all but ceased, we found that 6^/(1 + w) — > const, 
and that S^/{l + w) - 0{GM_2lR-2), where R-2 is the 
smallest radius for which 5 decreases as quickly as l/R^, 
M_2 is the mass inside R-2- 

Whilst it is possible with highly non-linear potentials, 
in general it seems improbably that 0(10^*^(1 4- w)) val- 
ues of which in turn result in order O(10~^(l -I- w)) 
values for will result in any pronounced non-linear 
effects especially over sub horizon scales. It is fair to 
say though that the conditions might be on the brink of 
breaking down for a 300/i^^ Mpc void, since HR ~ 0.1. 
But for smaller objects provided the conditions we place 
on the potential hold (which essentially preclude small 



changes in (j) from resulting in huge changes in V , V or 
V") then there is no particular reason to think that per- 
turbations in w will result in 0(1) or greater corrections 
to 5^. One should also point out that in the case where 
the scalar field is strongly coupled to niatter ( e.g. as 
in (iMota fc Shawl [20061 : iBrax et al.ll2nf)l iMota fc Sha^ 
12006D), then the dark energy perturbations might not 
be small. Since it is natural to expect that when mat- 
ter starts to infall into the nonlinear regime it will drag 
along the dark energy field. 

We also used our analytical expression for 5fi,/{l+w) to 
calculate the Integrated Sachs Wolfe (ISW) induced by 
the DDE density perturbations. We found that the ISW 
due to the DDE perturbations is always much smaller 
than the ISW effect due to the matter perturbation in 
these DDE models by a factor of roughly H^R^_2- Indeed, 
even the ISW effect due to the clustering of DDE near an 
extremely large 300 Mpc is found to result in a temper- 
ature anisotropy no larger than O(10~^(l + w)< 10^*). 
For smaller objects the ISW effect due to DDE clustering 
scales as (where Rc is the radius of the cluster). It 
is therefore unlikely that, if dark energy is described by 
a simple uncoupled quintessence theory, DDE perturba- 
tions will be detected through the ISW effect. 

In summary: we have derived an analytical expression 
for the DDE density contrast in uncoupled quintessence 
models and used it to make quantitative predictions for 
DDE clustering. Our results should also apply to fields 
with very weak matter coupling. This analysis could also 
be straightforwardly extended to models that are coupled 
to matter but we leave this to a later work. If the mass 
of the scalar in these models is 0{H), then 6^/(1 + w) 
is, to leading order, independent of the specifics of the 
underlying theory. As opposed to what has been sug- 
gested elsewhere, we found that 6^/(1 + w) is always 
small compared to the matter density contrast, and 
that even when the matter perturbation goes non-linear, 
S^/{1 + w) remains < 10~^ for typical astrophysical ob- 
jects. If we define the radius, Rc, of a matter over/under 
- density by the condition that at R = Rc-, the mean 
matter density contrast drops off like l/R^, then we 
found that \S^/d\ ~ 0((1 + w)H^Rl). DDE clustering 
could potentially be detected through the ISW effect or 
measurements of the peculiar velocity field of matter. 
Roughly, the magnitude of the DDE clustering contri- 
bution to these effects is a factor |(1 — ilm)Sci,/^mSm\ ~ 
0{{l-Q.„i){l+w)H'^ Rl/Q,rn) < 1 times the contribution 
from the clustering of ordinary matter. For realistic as- 
trophysical objects, we found that voids and super-voids 
correspond to local overdensities in dark energy, with 
S^/{1 + w) ^ 0(10^^) for clusters and average voids, 
and S^/{1 + w) ^ 0(10"^) for super-voids. If voids with 
radii of 100 — 300 Mpc exist within the visible Universe 
then Sij, may be as large as 10~'^(1 -I- w). Linear overden- 
sities of matter and super-clusters generally correspond 
to local voids in dark energy; for a typical super-cluster: 

V(i + ^) ~o(-io-')-^ 

Astronomical observations indicate w = — 1 ± 0.1 for 
2 < 1, and so we conclude that if dark energy is described 
by an uncoupled quintessence model, then today dark en- 
ergy is almost homogeneous on sub-horizon scales with 
perturbations generally of order of 10~^; the largest per- 
turbations will be 0(10^^ — 10"'^) and associated with 
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very large voids of matter. 
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